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Rosenfeld [Phys. Rev. A 15, 2545 (1977)] originally noticed that casting the transport coefficients 
of simple monatomic, equilibrium fluids in a specific dimensionless form makes them approximately 
single-valued functions of excess entropy. This observation has predictive value because, while the 
transport coefficients of dense fluids can be difficult to estimate from first principles, the excess 
entropy can often be accurately predicted from liquid-state theory. In this work, we use molecular 
simulations to investigate whether Rosenfeld's observation is a special case of a more general scaling 
law relating the tracer diffusivities of particles in mixtures to the excess entropy Specifically, we 
study the tracer diffusivities, static structure, and thermodynamic properties of a variety of one- 
and two-component model fluid systems with either additive or non-additive interactions of the 
hard-sphere or Gaussian-core form. The results of the simulations demonstrate that the effects 
of mixture concentration and composition, particle-size asymmetry and additivity, and strength of 
the interparticle interactions in these fluids are consistent with an empirical scaling law relating 
the excess entropy to a new dimensionless (generalized Rosenfeld) form of tracer diffusivity, which 
we introduce here. The dimensionless form of the tracer diffusivity follows from knowledge of the 
intermolecular potential and the transport/thermodynamic behavior of fluids in the dilute limit. The 
generalized Rosenfeld scaling requires less information, and provides more accurate predictions, than 
either Enskog theory or scalings based on the pair-correlation contribution to the excess entropy. As 
we show, however, it also suffers from some limitations, especially for systems that exhibit significant 
decoupling of individual component tracer diffusivities. 



I. INTRODUCTION 

Many computational and experimental studies have 
now provided empirical evidence of a strong correlation 
between transport coefficients and the excess entropy of 
equilibrium fluids (see, e.g., [3, i, i, H i, @, 0). The 
transport coefficients (e.g., diffusivity, viscosity, and ther- 
mal conductivity) quantify the dynamic response of a 
fluid to a small perturbation in the associated field vari- 
ables, while the excess entropy is a negative quantity that 
characterizes the number of microstates rendered inacces- 
sible to the fluid (relative to an ideal gas) due to static 
interparticle correlations. Changes to macrostate vari- 
ables that strengthen the interparticle correlations, and 
hence make excess entropy more negative, typically result 
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The connection between transport coefficients and ex- 
cess entropy is of fundamental interest because it pro- 
vides a clue in the long-standing puzzle concerning what 
structural and thermodynamic properties correlate with 
the dynamics of equilibrium fluids. The link also has 
practical consequences. For example, if the transport co- 
efficients of a fluid, cast in an appropriately reduced form, 
can be approximately represented as a single- valued func- 
tion of the excess entropy, then knowled ge o f the latter 
allows indirect "prediction" of the former [Tl| ■ The value 
of this approach lies in the fact that, while it is difficult to 
directly estimate transport coefficients from first princi- 
ples, the excess entropy can often be accurately predicted 
from liquid-state theories. 

At present, a rigorous and general statistical mechani- 
cal justification for the empirically observed relationship 
between transport coefficients and excess entropy is lack- 
ing. However, even in the absence of a formal justifi- 
cation, there are a number of practical questions that 
deserve further investigation. Here, we present new cal- 
culations that address key aspects of two such questions. 

• For what types of fluid systems is tracer diffusivity, 
when cast in an appropriately reduced form, ap- 
proximately a single-valued function of excess en- 
tropy? 

• Can we develop a strategy for determining the 
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aforementioned "appropriately reduced form" for 
tracer diffusivity of a given system from knowledge 
of the intermolecular potential, temperature, and 
composition? 

To understand the context of these questions, it is 
helpful to first consider some background. It has long 
been appreciated that the following reduced form of self- 
diffusivity, D R = Dp 1 / z {m/k-B,T) 1 / 2 ^ can be formally rep- 
resented as a single-valued function of excess entropy, 
s ex , for any model fluid of identical particles which in- 
teract via an inverse-power-law (IPL) pair potential of 
the form v(r) — e(a/r) n (see, e.g., [33j). Here, D is 
self-diffusivity, p is number density, T is temperature, fee 
is the Boltzmann constant, m is the particle mass, and 
the combination eer" is the single parameter of the IPL 
potential. The function relating Z? R and s cx strictly de- 
pends on the exponent n, but the dependence is weak. 
In fact, Rosenfeld first pointed out that the relation- 
ship is "quasi-universal" in the sense that, for a given 
value of excess entropy, there is less than 30% variation 
in the predicted self-diffusivitics for different equilibrium 
IPL fluids with 4 < n < oo. Because of this observa- 
tion, we refer to D R in this work as the Rosenfeld form 
of reduced diffusivity. Recently, Mittal et al. [i[ demon- 
strated that the same quasi-universal relationship also 
adequately describes the correlation between D R and s ex 
for both bulk (isotropic) and confined (inhomogeneous) 
equilibrium Lennard- Jones fluids. An important point 
of the Rosenfeld Q and Mittal et al. [jj investigations, 
discussed from a different angle in a more recent study 
by Dyre and co-workers [34| , is that since both the static 
and dynamic properties of many dense, simple liquids are 
dominated by their repulsive interactions, they closely 
mimic the behaviors of "equivalent" IPL fluids. Thus, 
one can expect D R of these monatomic simple-liquid sys- 
tems to scale with s cx in a way that is consistent with 
the trend originally identified by Rosenfeld. 

Yet there are many fluids that cannot be expected to 
mimic the behavior of monatomic IPL systems. Will 
excess entropy also prove useful for predicting the dy- 
namics of these more complex fluids? For example, can 
excess entropy be used to reliably predict the effects of 
temperature and density on the self-diffusivity of model 
fluids with soft (or even bounded) pair potentials, like 
those that characterize the effective interactions between 
macromolecules or micelles in solution [H, |H, At 
present the answer is unclear. A recent molecular dy- 
namics simulation study by Krekelberg et al. [3(| demon- 
strated that D R of the Gaussian-core fluid is not (even 
approximately) a single-valued function of s ex , and the 
same is true for the Rosenfeld-scaled thermal conduc- 
tivity and viscosity [3q . Is there a systematic way to 
construct, based on knowledge of the pair potentials and 
temperature of a system, an alternative reduced form, 
i.e. a generalized Rosenfeld diffusivity D GR , that (to 
within acceptable tolerances) is a function of excess en- 
tropy alone? One of the goals of this study is to address 
this question for different types of model systems with a 



variety of interactions, [e.g., continuous or discontinuous, 
steeply repulsive (diverging) or soft (bounded)] . 

A related question is whether excess entropy can be 
used to predict the effects of temperature, density, and 
composition on the tracer diffusivities of the components 
of a fluid mixture? This question has been recently stud- 
ied in a limited context. Specifically, following initial 
work on monatomic systems by Dzugutov a scaling 
for the tracer diffusivities of mixtures based on two-body 
contributions to the excess entropy has been introduced 
[39l | . Although it appears that this mixture version of 
the Dzugutov scaling can capture some of the behaviors 
exhibited by simple fluid systems, it also has some sig- 
nificant limitations. For example, the reduced diffusivity 
for the Dzugutov scaling, D z , relies on defining an effec- 
tive "hard-core diameter" for each interparticle potential, 
which is not convenient for the study of soft or penetra- 
ble particles with bounded interactions. Moreover, the 
Dzugutov scaling fails to describe the behavior of systems 
in the limit of vanishing number density. Finally, com- 
puting the two-body excess entropy requires knowledge of 
the radial distribution functions between all components 
in the mixture for each thermodynamic state of interest, 
which is particularly cumbersome when studying inho- 
mogeneous fluids. This should be contrasted with the 
excess entropy used in the Rosenfeld scaling, which can 
be readily calculated from knowledge of the fluid's equa- 
tion of state. For all of these reasons, we examine in this 
study whether one can, based on knowledge of mixture 
composition, temperature, and pair potentials, construct 
a generalized Rosenfeld form for the reduced tracer dif- 
fusivity for component i, -Dp R , that is approximately a 
single- valued function of the excess entropy of the fluid 
mixture. 

The organization of the paper is as follows. In Sec- 
tion nn we introduce the simple idea that underlies the 
generalized Rosenfeld scaling for predicting tracer dif- 
fusivity from excess entropy. Section [Ml provides details 
on the model fluid systems and the simulation techniques 
that we use here to put the predictions of this scaling to 
quantitative tests. In Section IIV1 we analyze the gen- 
eralized Rosenfeld scaling for a wide variety of binary 
mixtures of hard spheres, Widom-Rowlinson particles, 
and Gaussian-core particles. In Section [Vj we discuss 
how this data helps to understand the strengths and lim- 
itations of using excess entropy for predicting the effects 
that macroscopic parameters (temperature, density, com- 
position) and microscopic details (particle diameter, par- 
ticle mass, softness of the interparticle potential) have on 
single-particle dynamics. 



II. GENERALIZED ROSENFELD FORM FOR 
REDUCED TRACER DIFFUSIVITY 

In order to ensure that Dp R , the generalized Rosenfeld 
form of the reduced tracer diffusivity of component i, is 
defined in a way that is useful for making excess-entropy 
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based predictions, we aim to have it satisfy the following 
three criteria for a given system, (i) It should be propor- 
tional to the bare tracer diffusivity, i.e., -Dp R = <xDj. (ii) 
The prefactor a, which has units of reciprocal diffusivity, 
should be readily calculable based on knowledge of the 
parameters that define the fluid system, i.e., macroscopic 
variables like temperature T, density p, and the mole 
fractions of the species, as well as microscopic parame- 
ters like particle masses and the pair potentials Vij (r) de- 
scribing the effective interactions, (iii) The dimensionlcss 
quantity Df K should be approximately a single-valued 
function of s ex . 

For the case of a monatomic IPL fluid, the aforemen- 
tioned criteria are rigorously satisfied at all state points 
if one adopts the Rosenfeld form of the reduced self- 
diffusivity D R discussed in the introduction. In fact, 
Rosenfeld previously illustrated that, at low num- 
ber density, an analytical relationship between D R and 
s cx can be obtained for an IPL fluid by using two equa- 
tions: an Enskog theory expression for D and a truncated 
second- virial expansion for s cx . 

Our approach here is to similarly examine the low- 
density limit of more complex model fluids and mix- 
tures, taking advantage of known theoretical results for 
the tracer diffusivity and excess entropy to seek out a 
potentially useful definition for D^. As we show below, 
these low-density theoretical results do suggest a simple 
expression for Df K that, in the dilute limit, satisfies the 
three criteria mentioned above. Of course, unlike for an 
IPL fluid, choosing Z?p R of a complex fluid so that it is 
a single-valued function of s cx at low density does not 
guarantee that it will also behave that way at high parti- 
cle density. In fact, one can view the extent to which the 
Dp R versus s ex scaling holds at higher particle density as 
a measure of the utility of excess entropy for predicting 
tracer diffusivity of a given fluid system. Section [IV] fo- 
cuses on quantitatively examining this point for a variety 
of model fluids. 

To make the above discussion more concrete, first con- 
sider that kinetic theory indicates that Di is inversely 
proportional to number density p in a fluid mixture at 
low density [4(| 



Di = 



[pDi]o 
p 



(1) 



where the quantity [p-Di]o = lim p ^o pDi generally de- 
pends on temperature, mixture composition, as well as 
the masses of the species and their the interparticle in- 
teractions. We discuss simple theoretical methods for 
estimating [pDi]o for model systems below. The rela- 
tionship between p and s cx , to leading order in p, can be 
expressed as 
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dT 
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where B is the second virial coefficient, given by 



(3a) 



the sums are over the N c components of the mix- 
Xi is the mole fraction of component i, and By can 
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ture, 

be expressed in terms of the pair potentials Vij (r) as 
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where 1 = k^T. Using Eq. to eliminate p from 
Eq. |T]) and rearranging leads to 



Di 
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[pDi] -s cx /fc B 
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which again is valid only in the the p — » limit. We 
identify the dimensionlcss quantity on the left-hand-side, 
which is clearly a function of s ox only at low density, 
as the generalized Rosenfeld reduced form of the tracer 
diffusivity, L>f R , 



D 



GR 



Di 



(B + T^) [p Dl 



(5) 



Note that the expression for Z? GR , the generalized Rosen- 
feld self-diffusion coefficient for a monatomic fluid, is ob- 
tained by replacing Di with D in Eq. ([5]). 

Although the definition for Dp R given in Eq. ([5]) is 
compact, it is more convenient for making predictions 
for model systems if [p-Djo is further expressed in terms 
of the mole fractions of the species, the associated pair 
potentials, the particle masses, and the temperature. Be- 
low, we present simple theoretical expressions for carry- 
ing this out for particles with hard-sphere and soft (con- 
tinuous) interactions, respectively. 



A. Hard-Particle Interactions 

For models with hard-sphere interactions, an expres- 
sion for \pDi\p is easily obtained within Enskog kinetic 
theory [4l|,|42| In particular, the product pDi is given by 



pD t = 



87T 1 / 2 



-1/2 • 



(6) 



where rrn is the mass of component i, aij is the hard- 
sphere contact diameter between particles of type i and j, 
and gij (cr^ ) is the value of the radial distribution function 
between particles of type i and j at contact. The low 
density limit [p-E>i]o is obtained by substituting gij{cr^) — 
1 into Eq. ([6]), which gives 

[pDi]o = 



8vri/2 „ 

2^3=1 



2 v mj 



-1/2- 



(7) 
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When computing Z?p R for the hard-sphere and Widom- 
Rowlinson model mixtures discussed in Section IIIH we 
simply substitute Eq. (0 into Eq. ([5|). 



In Eq. ([10]) . Vij is the pair potential between particles of 
species i and j, and terms involving three-body static cor- 
relations have been omitted. In order to evaluate [pDi] Q , 
we take the p — » limit of Eqs. (JSJ)— (JTUJ) , which leads to 



B. Soft-Particle Interactions 

We also study fluids of soft particles in this work, 
i.e., particles with continuous and bounded interactions 
that cannot be treated as hard spheres with an effec- 
tive temperature-dependent diameter. In order to pre- 
dict [pDi]o for these models, we use an approximate the- 
ory due to Tankeshwar and co-workers O |44| , which we 
refer to as the Tankeshwar diffusion model (TDM). We 
have found that this basic theoretical approach strikes 
a reasonable balance between simplicity and accuracy. 
It has been shown to semi-quantitatively describe how 
temperature, composition, and density affect the tracer 
diffusivity of Lennard-Jones fluids, the one-component 
plasma, and Yukawa fluids [45j • We have also found that 
it approximately captures how temperature and density 
affect the diffusion coefficient of the Gaussian-core fluid. 

The details of the TDM are discussed extensively else- 
where [43|, 0, [45|. In short, it is based on an approxi- 
mate expression for the velocity autocorrelation function, 
and hence the tracer diffusivity via the Green-Kubo 
relation, for each component i of the fluid in terms of 
two parameters: the Einstein frequency u>i and a "jump- 
ing" frequency Tj. The values of these parameters are 
obtained by ensuring that the velocity autocorrelation 
functions satisfy some exact microscopic sum rules. 

Within this model, the product pDi is given by 
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of which follow by replacing gij(r) in Eq. (|10[) . with the 
Boltzmann factor of the pair potential, 



limay(r) = exv[-j3Vij{r)}. 



(12) 



When computing I?p R for the Gaussian-core mixtures 
discussed in Section IIII1 we substitute Eq. pip into 
Eq. ©. 



III. MODEL SYSTEMS AND SIMULATION 
METHODS 

As discussed in Section [H a key aim of this study is to 
investigate whether it is possible to construct an excess- 
entropy based strategy for predicting tracer diffusivity 
generic enough to be successfully applied to fluid mix- 
tures with either hard (impenetrable) or soft (penetra- 
ble) interparticle interactions. For our model systems, we 
choose familiar representations for both: the hard-sphere 
(HS) pair potential for the former and the Gaussian-core 
pair potential for the latter. 

The HS pair potential is discontinuous and athermal, 
assigning infinite energy to configurations that have par- 
ticle overlaps and zero energy to all others. It is thus 
represented as 



r < a 
r > a, 



y - 



(13) 



Here, cr^ is the contact diameter between particles of 
type i and j. We investigate several binary HS mix- 
tures with additive diameters, i.e., aij = (<7j + <7j)/2. In 
particular, we first study compositional effects on tracer 
diffusivity and excess entropy using a system composed 
of equimass (mi/rri2 — 1) particles with diameter ratio 
Ci/cr2 = 1-3. For this system, we examine mole frac- 
tions of component one in the range 0.1 < xi < 0.9. 
We also investigate the effects of diameter ratio by con- 
sidering particles with (J\j<j2 = 1.3,2.0,3.0,5.0. These 
latter studies are carried out at fixed composition x\ = 
0.1. All of the above systems are studied across a wide 
range of packing fractions ip = n(xial + xio\)j§ span- 
ning between 0.05 (dilute gas) and 0.5 (near the single- 
component HS freezing transition). 

We also consider a highly non-additive version of the 
binary HS mixture: the Widom-Rowlinson (WR) model 
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[461] . In this system, the contact diameter between par- 
ticles of the same type is zero (<7u = er 2 2 = 0), but the 
cross-diameter is finite a 12 — °2i = a and mi = m 2 . 
As might be imagined, this system exhibits entropically 
driven phase separation at sufficiently high density, which 
we avoid here by studying < pa 3 < 0.7. Since the 
model is symmetric with respect to the interactions, we 
can deduce global behavior by studying mole fractions in 
the range < x\ < 0.5. 

Finally, we study fluids composed of soft particles that 
interact via the bounded Gaussian-core pair potential 
[47| . given by 

Vg c (r) = e ij exp[-(r/a ij )% (14) 

where e.y and Oij are parameters that characterize the 
energy and length scale, respectively, of the interaction 
between particles of type i and j. For the simulations in 
this work, we truncate the interparticle interactions at a 
separation of 3.2c7y. We examine both single-component 
and two-component Gaussian-core fluids. For the latter, 
we adopt the same parameters used in a previous inves- 
tigation of the static structure and thermodynamics of 
that system (48|. Specifically, we assign 022 = 0.665o"n 

and = V {pxi + °22)A £ n = £ 22 and e 12 = 0.944eu 
(which encourages mixing), and mi = m,2- We inves- 
tigate the binary Gaussian-core fluid at compositions 
0.1 < xx < 0.9, temperatures 0.05 < k B T/e n < 0.4, and 
densities 0.05 < pu 3 i < 1.0. At some temperatures, the 
maximum density in this range is not an isotropic fluid 
due to propensity of the system to phase separate. We 
excluded from our analysis any state points that showed 
thermodynamic or structural indications of phase sepa- 
ration. 

To explore the dynamic properties of the above sys- 
tems, we perform molecular dynamics (MD) simulations. 
For the HS and WR mixtures, we use a standard event- 
driven algorithm (4(| . For the binary Gaussian-core fluid, 
the equations of motion are integrated using t he velocity- 
Verlet method [50| with time step St — O.Oby/mia^/en. 
All MD simulations are carried out in the microcanon- 
ical ensemble with N — 3000-5000 particle using a 
periodically-replicated cubic simulation cell with volume 
V, chosen in accord with the desired fluid density. Tracer 
diffusion coefficients D t are calculated by fitting the long 
time average mean squared displacement of the i type 
particles (Srf) to the Einstein relation 6Dit — (Srf). 
Note that, for the case of the monatomic fluid, this def- 
inition of the tracer diffusivity reduces to the self diffu- 
sivity. We perform multiple independent simulations at 
several state points for each model, and we find the rela- 
tive standard error in tracer diffusivities to be less than 
1%. 

Thermodynamic properties of the Gaussian-core fluid 
mixtures are computed using grand-canonical transition- 
matrix Monte Carlo (GC-TMMC) simulations. These 
simulations are conceptually equivalent to a series of 
semigrand simulations performed over a range of fluid 
densities stitched together using ghost insertion/deletion 



moves. Details of this method can be found elsewhere 
[5ll I52I ]. These simulations require fixed values of the 
activity, {£i,£;2}, volume, V, and temperature, T, as in- 
puts. The activity is defined as & = A^ 3 exp{p,i / k B T) , 
where pi and Aj are the chemical potential and the ther- 
mal de Broglie wavelength of component i, respectively. 
All GC-TMMC simulations for the Gaussian-core fluid 
mixtures reported here use a system volume of V = 343. 
For the activities, the values of In £j that we use span from 
32.63 at the lowest temperature to 12.24 at the highest 
temperature investigated. Thermodynamic properties at 
other values of activity are obtained via the histogram 
re-weighting technique. The primary quantity obtained 
from GC-TMMC is the particle number probability dis- 
tribution H(N; {pii}, V, T). From this, excess entropies 
are trivially calculated (see 10]). System size effects in 
excess entropy for the Gaussian-core mixtures are found 
to be negligle by comparing results to a scries of simula- 
tions using a smaller volume of V — 216. We confirmed 
that the equation of state (and hence excess entropy) 
of Gaussian-core fluid mixtures produced from the GC- 
TMMC simulations is statistically indistinguishable with 
that produced from molecular dynamics simulations. 

GC-TMMC calculations for the WR mixtures are per- 
formed using a system volume of V = 343. For this 
fluid, simulations are completed with activity values of 
£1 = £2 = 1 j an d histogram re- weighting is applied to ob- 
tain thermodynamic quantities at other values of activity. 
System size effects are examined by performing simula- 
tions over a limited density and composition range with a 
volume of V = 1000, and are also found to be negligible. 

The excess entropy data we present for the binary 
HS mixtures is calculated from the accurate Boublik- 
Mansoori-Carnahan-Starling-Leland (BMCSL) equation 
of state [H, H3]. As a check, we compared the BM- 
CSL values for compressibility factor and excess entropy 
against those obtained via molecular dynamics simula- 
tions for selected state points as a function of particle 
diameter ratio and packing fraction, and we found the 
agreement to be excellent. 

Finally, we also compare the results for the binary sys- 
tems to corresponding single-component systems. For 
the single-component HS system we use the data of [ll[ 
and for the single-component Gaussian-core fluid we use 
the data of [301 ] . 

IV. RESULTS AND DISCUSSION 
A. Hard-sphere mixtures 

1. Compositional effects 

We begin by investigating the effects of composition 
on mixtures of HS particles with size ratio £7i/f7 2 = L3. 
Figure [Ta| displays the tracer diffusion coefficients Di of 
the two components as a function of total packing frac- 
tion ip for several different compositions, indicated by 
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FIG. 1: (color online) Properties of the binary HS mixture with particle diameter ratio (cri/ai) = 1.3, equal mass, and a 
variety of compositions, (a) Tracer diffusion coefficients Di and (b) (negative) excess entropy — s ox as a function of packing 
fraction <p. (c) Rosenfeld Df and (d) generalized Rosenfeld -Dp R tracer diffusivities as a function of — s cx . Filled and open 
symbols denote component 1 (large) and 2 (small), respectively. The color of symbols denotes the mole fraction of component 
1, xi, specified in the legend of (a). The solid line in each figure is the result for the single-component HS system. The dashed 
red line in (d) represents a least-squares fit of the data to Eq. (|15[) . which results in a = 0.95, A = 1.85, and /3 = 0.74. In (c) 
and (d), the insets are the same as the main plots, but on a log-log scale. 



the mole fraction of large particles, x\. As must be the 
case, when one of the species is present in high concen- 
tration, its tracer diffusivity approaches the value of the 
self-diffusion coefficient D of the single-component HS 
fluid at the same packing fraction if. However, the tracer 
diffusivity of the dilute component is generally different 
than D. In particular, when component 1 (the larger par- 
ticles) is dilute, one should expect D\ < D. This logic 
can be qualitatively rationalized by the fact that, on aver- 
age, motion of the larger solute would require larger local 
structural rearrangements (i.e., fluctuations) than for the 
motion of the smaller solvent particles. Conversely, by an 
analogous argument, one expects D 2 > D if component 
2 (the smaller particles) is dilute. The data in Figure [Tal 
is consistent with these expectations. 

It is interesting to note that the compositional varia- 
tion in lnZ?i is fairly insensitive to the value of (p. More- 
over, the excess entropy [Figure lib] exhibits almost no 
compositional dependence whatsoever, and its packing 
fraction dependence for any particular composition is 
nearly identical to that of the single-component HS fluid. 
All of this suggests that an appropriate composition de- 
pendent rescaling of the tracer diffusivity data might (ap- 
proximately) make it a single-valued function of excess 



entropy. Indeed, Figures [Tc] and [Id] show that the gen- 
eralized Rosenfeld tracer diffusivities -Dp R of Eq. [5] for 
both species collapse onto a single curve [that describing 
the single-component (SC) fluid data, Dgc \ when plot- 
ted versus excess entropy, while no data collapse occurs 
if the tracer diffusivities are naively represented in the 
original Rosenfeld reduced form Df . 



The single-component relationship D 



SC 



can be 

described by a piece-wise function. The form of its low- 
density (low — s GX ) scaling is an inverse power law, given 
by Equation ^ . From Figure lldl we infer that when 
—s ex /kB ^ 1, the relationship becomes approximately 
exponential. A least squares fit assuming these generic 
function forms, i.e., 



a[-s cx //j B ] _1 
Aexp[-Bs cx /k B ] 



-s ex /fc B < 1, 
-s sx /k B > 1, 



(15) 



yields a — 0.95, A = 1.85, and (3 = 0.74, and describes 
the simulation data very well (see red dashed line in 
Fig. [Id]). 

The data collapse of Figure [Tdl suggests that tracer dif- 
fusivities, D± and D 2 , of this mixture might also be ade- 
quately predicted using Eq. (|15|) together with knowledge 
of the pair potentials, composition, and excess entropy 
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FIG. 2: (color online) Ratio of tracer diffusivity predicted 
from (a) Enskog theory [Eq. (|6}] and (b) from excess entropy 
and the single-component HS result (generalized Rosenfeld 
scaling) [Eq. (I16|l ]. for a HS mixture with (oi/<72) = 1.3, 
equal mass, and a variety of compositions. Red dashed lines 
represent 20% relative error of prediction. Symbols have the 
same meaning as in Figure [1] 



of the mixture. Specifically, the generalized Rosenfeld 
scaling prediction for tracer diffusivity of component i is 
given by 



D 



predicted 



( s cx ) = [ P A]o 



B + T- 



JB 

It 



D. 



GR/„cxn 

sc V s ) 



(16) 



with B from Eq. ©, [pDi} from Eq. ©, and Df§-(s ex ) 



from the fit of the single-component data to Eq. ([15]) . 

One way to quantitatively assess the relative predictive 
ability of Eq. (fl6|) is to compare it to the results of, e.g., 
the basic Enskog theory given by Eq. J^J). Both equations 
require as inputs several pieces of information, including 
the form of the pair potentials and the mixture composi- 
tion. While predictions based on the generalized Rosen- 
feld scaling also require knowledge of s cx for the mixture 
and properties of the single-component system, Enskog 
theory requires knowledge of the state dependent contact 
values of the three partial radial distribution functions of 
the mixture. 

Figure [2] displays the relative error in the pre- 
dicted tracer diffusivity to the simulated diffusivity 
Df redlctcd / Di . Enskog theory [Figure [2a] provides good 
predictions at high values of tracer diffusivity (i.e., low 
ip). However, as the value of Di decreases (i.e. , <p 
increases) Enskog theory first underpredicts, then ulti- 
mately significantly overpredicts When looking at 
the entire range of <p studied here, 80% of the Enskog 
theory predictions lie within 20% of the molecular simu- 
lation data. On the other hand, the excess entropy based 
expression of Eq. (fT6|) [Figure [2b] predicts the tracer dif- 
fusivities semi-quantitatively for all state points investi- 
gated here (100% of predictions within 20% of the simu- 
lation data). 



2. Particle-size asymmetry effects 



We also study the effects of particle-size asymmetry 
on the relationship between excess entropy and tracer 
diffusivity by examining a scries of binary HS mixtures 
at composition x\ — 0.1 and packing fractions in the 
range < ip < 0.5. Particles of types 1 and 2 were 
taken to have identical masses, but we investigated sev- 
eral systems with different diameter ratios \a\jcri = 1.3, 
2.0, 3.0 and 5.0]. Figure [3al displays the tracer diffu- 
sivities, D\ and D2, for these systems. Increasing the 
magnitude of the diameter ratio leads to progressively 
larger deviation of the tracer diffusivities from the self- 
diffusion coefficient D of the HS fluid at the same ip. As 
expected, larger particles diffuse slower than smaller par- 
ticles [D\ < D < D^\. How is the excess entropy affected 
by increasing the ratio of particle diameters? Figure I3bl 
shows that increasing axfui at fixed tp and x\ systemat- 
ically decreases — s ox (i.e., weakens the static interparti- 
cle correlations) . This effect is qualitatively connected to 
the more efficient packing arrangements that spheres can 
sample when significant polydispersity is present [H, IHl • 

Given that increasing size ratio uniformly reduces 
structural correlations, but impacts the dynamics of large 
and small particles in different ways, it might not be 
surprising that tracer diffusivity data represented in the 
original Rosenfeld form, Df-, does not collapse when plot- 
ted versus excess entropy [Fig. [3c]. The Df- data for 
the smallest size ratio a\j<y% = 1.3 is qualitatively sim- 
ilar to the single-component result. However, there is 
significant deviation for larger diameter ratios, with ex- 
cess entropy underpredicting the mobility of smaller par- 
ticles and overpredicting that of larger particles. Fig- 
ure [3d] shows, however, that tracer diffusivity reduced in 
the generalized Rosenfeld form, Dp R , mostly collapses 
to the single-component curve when plotted versus — s cx . 
The most pronounced deviations are for the largest size 
ratio (cri/cr2 = 3, 5) at the highest packing fractions 
(0 > 0.45). 

Figure[4]quantitatively compares the predictions of En- 
skog theory [Eq. ©] with those based on the generalized 
Rosenfeld scaling [Eq. HU)]. At high values of Di (low 
ip), both methods provide accurate predictions. As be- 
fore, for decreasing Di (increasing ip), Enskog theory first 
underpredicts and then ultimately overpredicts the tracer 
diffusivities. The excess entropy based predictions never 
underpredict, but eventually overpredict the mobility at 
high values of ip. As an overall measure, the Enskog and 
the excess entropy expressions predict 70% and 80% of 
the tracer diffusivities within 20% of the simulated val- 
ues, respectively. Moreover, we note that while the excess 
entropy method predicts the tracer diffusivities of the two 
components with similar reliability, the Enskog expres- 
sion does well for the small particles (90% within 20%), 
but poorly for the large particles (50% within 20%). 
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FIG. 3: (Color online) Properties of the binary HS mixture at composition xi = 0.1, equal mass, and several size ratios (ai/az) 
[see legend in (a)], (a) Tracer diffusion coefficient D; and (b) (negative) excess entropy — s cx versus packing fraction tp. (c) 
Tracer diffusivity reduced in the original Rosenfeld Df~ and (d) generalized Rosenfeld Dp H forms as a function of — s cx . Filled 
and open symbols represent large (component 1) and small (component 2) spheres. Symbol shapes denote different particle 
diameter ratios, as described in the legend of (a). The solid line in each figure is the result for the single-component HS system. 
In (c) and (d) the insets are the same as the main plots, but on a log-log scale. 



3. Two-body excess entropy scaling 

As noted in Section [H an alternative excess entropy 
based scaling for diffusion was introduced by Dzugutov 
Q , who found that an appropriately reduced form of the 
self-diffusion coefficient D z for atomic fluids at moder- 
ate densities is nearly a universal function of the two- 
body contribution to the excess entropy s^ 2 \ Subse- 
quently, others have suggested a generalization of this 
scaling (39| to predict tracer diffusivities of fluid mix- 
tures. In the generalization, the reduced tracer-diffusion 
coefficient, defined as Df = Di/xi, where 



1/2 



is thought to approximately scale with the i-component 
contribution to the two-body excess entropy, defined 
as 39] 



,(2) 



1 N ° f 

/k B = --p^Xj / dr 

i=i J 



Note that the two-body excess entropy per particle is 

.(2) 



57] 



{ 5ij (r)ln 5ij (r)-[ 9 y(r)-l]}. (18) 



given by s (2) = J2i 

In Figure [5al and Ibl we examine the mixture generaliza- 
tion of the Dzugutov scaling for the HS systems discussed 
above. Specifically, we show data with fixed diameter ra- 
tio and varying composition in panel (a) and fixed compo- 
sition and varying diameter ratio in panel (b). Both sets 
of data more or less track the scaling. However, devia- 
tions from the single component curve appear systematic. 

(2) 

The single-component relation with s\ overpredict the 
small sphere mobility and underpredict the large sphere 
mobility. 

Unfortunately, the predictive value of this type of scal- 
ing is inherently limited by the fact that the single- 

(2) 

component data cannot access the large values of — 
realized by the large spheres in a mixture. The for- 
mer reach a value of 5.5&B at <p = 0.5, while the latter 
are greater than 10/cb for the largest size ratios exam- 
ined here. As a result, "predicting" tracer diffusivities 
of a mixture would require some systematic way of ex- 
trapolating the single-component curve by a substantial 
amount. As discussed in Section HI the scaling is also 
limited to systems, like the HS fluid, for which the inter- 
particle repulsions are steep enough to define an effective 
hard-core diameter to each interaction. Thus, it will be 
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FIG. 4: (color online) Ratio of tracer diffusivity predicted 
from (a) Enskog theory [Eq. (|6}] and (b) from excess entropy 
and the single-component HS data (i.e., generalized Rosen- 
feld scaling) [Eq. (|16f) ] . for a HS mixture with equal mass, 
composition x\ = 0.1, and a variety of size ratios {a\ja-z). 
Red dashed lines represent 20% relative error of prediction. 
Symbols have the same meaning as in Figure [3] 



of little use for studying systems with bounded interac- 
tions like the Gaussian-core potential or a other models 
that characterize the soft effective interactions between 
macro- or supramolecular species in solution [3^ |. 



B. Widom-Rowlinson mixtures 

Here we examine the behavior of the Widom- 
Rowlinson (WR) model fluid introduced in Section IIIII 
Recall that it is defined as a mixture of non-additive hard 
spheres with o\\ — 022 = but oyi — o. Figure [Sa| dis- 
plays the tracer diffusivity Di as a function of density for 
a several compositions x%. Note that the Di is always 
greater then the self-diffusion of the single-component 
HS fluid, since the number of collisions per unit time 
will clearly be less in the WR fluid than in the HS fluid 
at the same density. At x\ — 0.5, D\ = D2, since the 
fluid is symmetric. At fixed density, as x\ decreases, D\ 
decreases while D2 increases. This is because the dilute 
species will experience many more collisions per unit time 
(it has more neighbors of the opposite type) than the con- 
centrated species. Likewise, Figure [6b] shows — s ex for the 
WR fluids is always less than that of a single-component 
HS fluid of the same density. This is expected since par- 
ticles of the same type do not directly exclude volume 
from one another, which in turn reduces the entropic 
driving force for forming strong interparticle correlations. 
Decreasing x\ from 0.5 toward zero at fixed density de- 
creases — s GX because it increases the number of particles 
in the system that do not interact. 

As was the case for the HS fluid mixtures, Figure Irjcl 
shows that the tracer diffusivity reduced in the original 
Rosenfeld form, Df-, is not even approximately a single- 
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FIG. 5: (color online) Tracer diffusivity reduced in gener- 
alized Dzugutov form Df discussed in text versus (negative) 

(2) 

i-component contribution to two-body excess entropy — s\ of 
binary HS mixtures, (a) Particle diameter ratio o\jo2 = 1.3, 
equal mass, and a variety of compositions. Symbols have same 
meaning as Figure [JJ (b) Composition x\ — 0.1, equal mass, 
and a variety of size ratios. Symbols have same meaning as 
Figure [3] Insets are the same as the main plots, but with a 
log-log scale. 



valued function of s GX . On the other hand, Figure [6d] 
shows that the tracer diffusivity data cast in the general- 
ized Rosenfeld form, Df R , largely collapses when plotted 
versus excess entropy. Interestingly, Dp R of the WR fluid 
is well described by the mathematical form of the single- 
component HS data. The quality of the collapse is more 
easily seen on a log- log scale (inset to Figure l6d|) . 

Figure [7a| and [b] compare the accuracy of predicting 
tracer diffusivity of the components of the WR fluid 
based on Enskog theory [Eq. ©] versus excess entropy 
of the WR mixture and the single-component relation for 
the HS fluid [Eq. flB])]. As for the HS fluid, Enskog the- 
ory predicts 80% of the WR data within 20% of the sim- 
ulation values. In contrast, the excess entropy method 
predicts 97% of the data within 20% of the simulated 
tracer diffusivities. 

Lastly, since the WR model is composed of (non- 
additive) hard particles, it represents another good test 
case for the mixture generalization of the Dzugutov scal- 
ing. Figure [8l however, clearly shows that this two-body 
scaling does not collapse the WR data. In general, par- 
ticles of type i diffuse considerably faster than would be 
predicted based on the single-component HS fluid behav- 
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FIG. 6: (color online) Properties of the Widom-Rowlinson mixture, (a) Tracer diffusivities Di (i = 1, 2) and (b) (negative) 
excess entropy — s cx as function of density pa 3 , (c) Tracer diffusivity reduced in original Rosenfeld form Df^ and (d) generalized 
Rosenfeld form Dp R as a function of — s cx . Filled and open symbols denote component 1 and 2, respectively. The symbol type 
denotes the mole fraction of component 1, xi, indicated in the legend of (a). The solid line in each figure is the result for the 
single-component HS fluid. In (d) the insets are the same as the main plots, but on a log-log scale. 
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FIG. 7: (color online) Ratio of tracer diffusivity predicted 
from (a) Enskog theory [Eq. (|6}] and (b) from excess entropy 
and the single-component HS result (generalized Rosenfeld 
scaling) [Eq. (|16[l ], for the Widom-Rowlinson mixture. Red 
dashed lines indicate 20% relative error of prediction. Sym- 
bols have the same meaning as in Figure |S] 



ior and the i-component of the two-body excess entropy 
in the WR mixture. Moreover, the magnitude of the un- 
der prediction depends sensitively on composition. This 
breakdown of the mixture version of the Dzugutov re- 



lation for non-additive HS fluids indicates that it is not 
as widely applicable, even within the limited class of HS 
model fluids, as the generalized Rosenfeld scaling intro- 
duced here. 



1(T 




FIG. 8: (color online) Tracer diffusivity reduced in the mix- 
ture generalized Dzugutov form discussed in text Df versus 
the (negative) i-component of the two-body excess entropy 
sf 2 ' for the Widom-Rowlinson model. Symbols have same 
the meaning as Figure [6] 



11 



C. Single-component Gaussian-core fluid 

As discussed in Section HI and much more extensively 
in [33, [3ll, [HI, Hil , the properties of the single-component 
Gaussian-core fluid are anomalous compared to those of 
simple atomic liquids. For example, as shown in Fig- 
ure [9a] the self diffusivity D of the Gaussian-core fluid 
first decreases, and then anomalously increases as a func- 
tion of particle density along an isotherm. Likewise, — s cx 
at constant temperature first increases (the fluid becomes 
more structured), and then anomalously decreases (the 
fluid becomes less structured) as a function of density. 

In brief, these unusual trends can be qualitatively ra- 
tionalized based on the Gaussian form of the repulsion. 
When the density and temperature are sufficiently low, 
the distance between particles is larger than the range 
of the potential. Under these conditions, the part of the 
interaction that the particles sample when they "collide" 
appears steeply repulsive, and thus the effects of density 
on structure and dynamics are similar to those of HS 
fluid. However, at high particle densities, particles in the 
Gaussian-core fluid effectively overlap one another due to 
the bounded form of the interaction. The effect is that 
each particle constantly experiences largely canceling soft 
repulsive forces of many neighbors. Increasing the par- 
ticle density under these conditions enhances this effect, 
paradoxically weakening the structural correlations and 
increasing the self diffusivity of the fluid. 

It is clear from Figure [5a| and [b] that both D and 
s cx are strongly correlated for the Gaussian-core fluid. 
In Figure [Sc] we show the self-diffusion coefficient ex- 
pressed in the original Rosenfcld form, D R , as a function 
of — s cx . As noted previously [3(|, this basic scaling is 
not even approximately a single-valued function of ex- 
cess entropy. However, similar to the behavior of the HS 
and WR mixtures discussed previously, Figure [9d] shows 
that the generalized-Rosenfeld-scaled self diffusivity Df K 
collapses to a single curve when plotted versus excess en- 
tropy. The quality of the collapse at even low density is 
apparent in the log-log plot shown in the inset to Fig- 
ure I9dl We also observe that the same functional form 
that was used to fit the single-component HS data can 
also be applied to the Gaussian-core system. In particu- 
lar, fitting the data to Eq. ([15]) yields a = 0.59, A = 1.33, 
and B = 0.90. As shown in Figure Ed] Eq. ([15]) with 
these parameters (red dashed line) describes the simula- 
tion data very well. 



D. Binary Gaussian-core mixtures 

As a final test of the relationship between single- 
particle dynamics and excess entropy in soft-particle flu- 
ids, we examine the binary mixture of Gaussian-core par- 
ticles described in Section uTTl In particular, Figures [TOaT - 
[b] display the tracer diffusivities of the large and small 
Gaussian-core particles as a function of density for a va- 
riety of mixture compositions (x± — 0.1, 0.3, 0.5, 0.7, 



and 0.9) and reduced temperatures k^Tjt — 0.05, 0.1, 
and 0.2. The first point of interest in the data, evident in 
Figure fTUal is that the tracer diffusivity of the larger type 
1 particles displays the same anomalous trend as a func- 
tion of density as the single component Gaussian-core 
fluid. That is, increasing the density eventually leads to 
an anomalous increase in D\. However, over the density 
range considered here, the tracer diffusivity of the small 
type 2 particles does not show this anomalous trend (see 
Figure [10b}. From a qualitative perspective, these differ- 
ent behaviors perhaps might be expected, since the larger 
particles begin to overlap more (and hence transition into 
anomalous, mean-field behavior) at lower densities than 
the smaller particles. This aspect of binary Gaussian- 
core mixtures in discussed in detail elsewhere |3l| . 

One consequence of the dynamic decoupling of small 
and large particles described above is a crossover density 
for tracer diffusivity. Specifically, large particles have 
lower tracer diffusivity than small particles at low den- 
sity, but they attain higher values of tracer diffusivity 
than small particles at sufficiently high density [see Fig- 
ure [10c]. Because of this crossover, the Gaussian-core 
fluid mixture is an interesting counterexample to the flu- 
ids discussed thus far. It appears that both components 
cannot scale in a simple way with a single static measure 
like — s cx [Figure [lOdj . This is evident, when one con- 
siders how the reduced Rosenfeld [Figure UOdj and gen- 
eralized Rosenfeld [Figure HOej forms of tracer diffusiv- 
ity behave as a function of — s cx . As before, the original 
Rosenfeld form, Df, fails to collapse any of the data. The 
generalized Rosenfeld form, -Dp R , does an excellent job 
of collapsing the low-density data, but necessarily breaks 
down at higher densities, where the anomalous behavior 
emerges. 

We close our discussion of the binary Gaussian-core 
mixture with an interesting empirical observation. In 
the spirit of [oOf , we find that a collective tracer diffu- 
sivity of the mixture, which we define here as D GK C = 
(Df R ) Xl (D2 R ) X2 , is still a single- valued of excess entropy 
over the wide range of temperature, density, and compo- 
sitions investigated here [see Figure [TUg]. As can be seen, 
it also quantitatively tracks the relationship between D R 
and — s ox for the single-component Gaussian-core fluid. 
What this implies is that tracer diffusivity of one com- 
ponent can be predicted based on knowledge of tracer 
diffusivity of the other component, the excess entropy 
of the mixture, and the behavior of the one-component 
fluid. Of course, this observation also holds true (triv- 
ially) for the other mixtures we discussed earlier because 
the generalized Rosenfeld tracer diffusivities themselves 
are single- valued functions of s cx of those systems. 



V. CONCLUSIONS 

In this work, we present a new dimensionless form of 
the tracer diffusion coefficient of a species which we call 
the generalized Rosenfeld tracer diffusivity. We show, 
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FIG. 9: (Color online) Properties of the single-component Gaussian-core fluid, (a) Self diffusivity D and (b) (negative) excess 
entropy — s cx versus density p. (c) Self diffusivity reduced in the Rosenfeld form D R and (d) the generalized Rosenfeld form 
D Gn as a function of — s cx . In (a) and (b), arrows indicate increasing temperature T. In (c) and (d), symbol type corresponds 
to density, indicated in the legend of (c). In (d), inset is the same as the main plot but on log-log scale. In (d), the red dashed 
line represents a least-square fit of the data to Eq. (|15[1 . 



via molecular simulation, that this quantity is approxi- 
mately a single-valued function of excess entropy for a 
range of model one- and two-component fluid mixtures. 
The empirical excess entropy scaling is consistent with 
the various effects that composition, temperature, den- 
sity, and microscopic interactions have on the equilibrium 
single-particle dynamics of these systems. Generalizing 
an earlier argument of Rosenfeld [2j, we show that the 
functional form of the reduced tracer diffusivity can be 
obtained by examining the theoretical behavior of excess 
entropy and tracer diffusivity in the low-particle-density 
limit. 

We demonstrate that the aforementioned "generalized 
Rosenfeld" scaling applies more broadly than other sim- 
ple approaches such as Enksog theory or empirical scal- 
ings based on the pair-correlation contribution to the ex- 
cess entropy. However, we also identify some important 
limitations of the approach. For example, the scaling 
breaks down for highly asymmetric hard-sphere mixtures 
(diameter ratios of 5 or larger) for packing fractions near 
the freezing transition. It also breaks down for Gaussian- 
core mixtures, where the softness of the interactions com- 
bined with the size asymmetry gives rise to significant 
decoupling of the single-particle dynamics of the species. 
Interestingly, even in this latter case, we show that a sin- 
gle, collective measure of the tracer diffusivities obeys an 



excess entropy scaling, which provides a quantitative link 
between structure and the tracer diffusivities of the two 
components. 



It may also be fruitful in future work to develop gener- 
alized Rosenfeld scalings for other transport coefficients, 
like thermal conductivity and shear viscosity. We plan to 
focus on extending the ideas of the present study to sys- 
tems with other types of dynamics (e.g., including effects 
of dissipation, hydrodynamic interactions, etc.). We also 
are studying what aspects of interparticle interactions 
can give rise to decoupling of species-specific structural 
and dynamic quantities. 
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FIG. 10: (color online) Properties of the binary Gaussian-core fluid described in the text. Tracer-diffusion coefficient of (a) 
component 1 and (b) component 2 versus density, (c) Illustration of the crossover in tracer diffusivity for components 1 and 
2 as a function of density, (d) Excess-entropy — s ox versus density, (e) Rosenfeld scaled tracer-diffusivity Df^, (f) generalized 
Rosenfeld tracer diffusivity -D S GR , and "collective" generalized Rosenfeld tracer diffusivity D GR ' C = (D^ K ) X1 {D^Y 2 versus 
— s cx . Symbol type corresponds to reduced temperature hsT/en: 0.05 (circles), 0.1 (squares), and 0.2 (diamonds). For clarity 
in (a)-(d), increasing temperature is given by the direction of the arrow.. In (a)-(f), closed and open symbols denote components 
1 and 2, respectively. Panel (c) displays the crossover behavior of the tracer diffusivities for xi = 0.9 and ksT/en — 0.05 and 
0.1. In panels (f) and (g), small black crosses represent the single-component Gaussian-core data. Inset to (g) is the same as 
the main plot but on a log-log scale. 



APPENDIX: EXTENSION TO BROWNIAN 
DYNAMICS 

Consider a collection of Brownian particles of radius a 
and volume fraction ip = 47ra 3 p/3 suspended in a contin- 
uum solvent. We wish to identify a generalized Rosenfeld 
scaling of D, the long-time self-diffusion coefficient of the 
particles, which will be (approximately) a single-valued 
function of the excess entropy. In this case, since the 
solvent can be being treated as a continuum, the excess 
entropy of interest is that associated with the static cor- 
relations of the Brownian particles. 

At infinite dilution, the excess entropy of the Brownian 
particles is zero, and D is simply equal to the Stokes- 
Einstein diffusivity Dq = k^T / (Sirr]a) where 5 = 4 for 
slip and S = 6 for stick boundary conditions, and rj is the 
solvent viscosity. To leading order in packing fraction 
ip, the difference between Dq and D can be expressed 



[6H,|6a|63j,|6 



AD ee D - D = -D D 2V , 



(A.l) 



where Di characterizes how static correlations modify 
the long-time self-diffusivity of the particles. This quan- 
tity can be expressed as 

f-OO 

D 2 = {-Z + A 11 +2B 11 )g{r)r 2 dr+ 
Jo 

An — A12 — Bn + B12 



1 ( dAu 

2 \ dr 



dA 



12 



dr 



Q(r)g(r)r 2 dr. (A.2) 



Expressions for the functions Q(r), An, A^, Bn, and 
B12 are known in the hydrodynamic limit 62],[65j|. Also, 
in the dilute limit, one may further replace g(r) in 
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Eq. (|A~2|) by the Boltzmann factor [Eq. ([T2j)] of V'(r), 
the effective pair potential between the Brownian parti- 
cles in solution. 

Following Eq. one can similarly express the ex- 
cess entropy s' cx associated with the structure of the 
Brownian particles (to leading order in particle den- 
sity) in terms of its osmotic second virial coefficient, 
s' cx = -k B p[B' + dB'/dlnT]. As in Eq. (3b}, B' can 
be readily obtained from knowledge of V'(r). Using the 
osmotic virial expression for the excess entropy to elim- 
inate packing fraction from Eq. (jA.lj) and rearranging 



yields a Brownian-appropriate generalized Rosenfeld self 
diffusivity: 



AD GR = 



(D - D) 



B' + T% 



-DnD 



V cx /fc B . (A.3) 



Although the quantity, AD GK is strictly a single-valued 
function of s' cx in the dilute limit, its behavior at higher 
particle concentrations needs to be studied further. 
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